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Abstract 

We systematically study the possibility for realizing realistic values 
of lepton mass ratios and mixing angles by using only renornralizable 
Yukawa couplings derived from the heterotic Zq-1 orbifold. We as¬ 
sume one pair of up and down sector Higgs fields. We consider both 
the Dirac neutrino mass scenario and the seesaw scenario with degen¬ 
erate right-handed majorana neutrino masses. It is found that realis¬ 
tic values of the charged lepton mass ratios, m e /m T and m^/m T , the 
neutrino mass squared difference ratio, and the lepton 

mixing angles can be obtained in certain cases. 
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1 Introduction 


Understanding the origin of fermion masses and mixing angles is one of the 
most important issues in particle physics. Quark masses and mixing angles, 
and charged lepton masses are well known. They have hierarchical structures. 
Recently, neutrino oscillation experiments have provided us with information 
about neutrino mass squared differences and mixing angles [T] . Two of these 
mixing angles are large in contrast to the small quark mixing angles. Within 
the framework of the standard model and its extensions, a fermion acquires 
its mass from the Yukawa coupling with the electroweak Higgs fields. For 
right-handed neutrinos, majorana masses are additional sources of masses 
and mixing angles. Thus, it is important to study Yukawa couplings and 
right-handed majorana neutrino masses derived from an underlying theory. 

Superstring theory is a promising candidate for a unified theory including 
gravity. For consistency, superstring theory requires 10D space-time. That 
is, it predicts 6D extra space other than our 4D space-time, and such a 
space must be compact. This 6D compact space is very important from the 
viewpoint of string phenomenology. In general, it is an origin of fermion 
flavor structure. That is, the flavor structure derived from a string model 
depends on the geometrical aspects of the 6D compact space, and in principle 
one can calculate Yukawa couplings for a given compact space with known 
geometry. Geometrical aspects of the 6D compact space also provide us with 
selection rules for allowed Yukawa couplings. In several cases, such selection 
rules are so severe that off-diagonal Yukawa couplings are not allowed for one 
Higgs held, because different families are discriminated by quantum numbers, 
which originate from geometrical aspects of the 6D compact space. 

Among several types of string models, heterotic orbifold models J2] and 
intersecting D-brane models are particularly interesting, 4 to realize realistic 
Yukawa couplings. They have localized modes in the 6D compact space, and 
their Yukawa couplings can be calculated 0 0 ei m , because their 6D geom¬ 
etry is not complicated. 5 Indeed, a Yukawa coupling among localized modes 
has a suppression factor depending on their distances 0 0 HH El 0 • That 
can explain suppressed Yukawa couplings. Therefore, it is very important 
to study possibilities for deriving realistic fermion masses and mixing angles 
from heterotic orbifold models and intersecting D-brane models. The number 

4 See for a review of intersecting D-brane models, e.g. Ref. [3| and references therein. 

5 See for calculations of Yukawa couplings in intersecting D-brane models Refs. HOE! EDI- 
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of 6D Z N orbifolds (2, E21E3] and Z N x Z M orbifolds t T3J, especially with 4D 
N — 1 supersymmetry, is finite, i.e. Z 3 , Z 4 , Z 6 -I, Z 6 - II, Z 7 , Z 8 -I, Z 8 - II, Z 12 -I, 

Z 12 -II, and Z N xZ M for (. N , M) = (2, 2), (2, 3), (2,4), (2, 6), (3, 3), (3, 6), (4,4), (6,6). 
(See Refs, [^j for examples of explicit models, and Refs. ESI E3 EH EH] for 
recent model buildings.) In particular, non-prime order orbifolds would be 
interesting, because they allow off-diagonal couplings HU EH- On each orb- 
ifold, all the fixed points are known, where light modes can be localized. 

Thus, a systematical analysis on Yukawa matrices, which can be derived 
from heterotic orbifold models, would be possible, while the number of inter¬ 
secting D-brane configurations seems to be infinity. Indeed, such studies have 
been done for the quark sector [2211231 • In particular, in Ref. B3 possibili¬ 
ties for leading to the realistic mixing angle V c b and mass ratios m c /m t and 
m s /mb from the Zq-1 orbifold have been shown in the case with the minimal 
number of Higgs fields. 

However, such a study considering mixing angles has not been done for 
the lepton sector. In this paper, we study lepton mass ratios and mixing 
angles derived from the Z 6 -I orbifold. Commonly, the smallness of neutrino 
masses is explained in two ways. One is the Dirac neutrino mass scenario, 
that is, Yukawa couplings between neutrino and Higgs fields are strongly 
suppressed for some reason. Since a Yukawa coupling is suppressed depend¬ 
ing on the size of extra dimensional space in an heterotic orbifold model, 
strongly suppressed Yukawa couplings can be obtained in the case that an 
extra dimensional space is large compared with the string scale. The other 
is the seesaw scenario [23], which requires right-handed majorana neutrino 
masses at the intermediate scale between the weak scale and the Planck 
scale. 6 Within the framework of string theory, the natural mass scale is just 
the string scale, and such an intermediate mass scale should be obtained 
through vacuum expectation values (VEVs) of some scalar fields. However, 
that is quite model-dependent. Here we study both the Dirac scenario and 
the seesaw scenario, and for the latter case we assume the right-handed ma¬ 
jorana neutrino mass matrix to be proportional to the identity matrix with a 
universal mass scale. With this, we consider all of the possible assignments of 
leptons to fixed points, examine Yukawa terms allowed by the selection rule 
and calculate their Yukawa matrices varying orbifold moduli parameters. In 
practice, the number of relevant moduli parameters is two in our models. 

6 In Ref. (25) , realization of the minimal seesaw mechanism has been examined in explicit 
Z 3 orbifold models, and its difficulty has been shown. 
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Then we try to fit the Yukawa matrices by these two moduli parameters in 
order to get realistic values of six observables, that is, lepton mass ratios and 
mixing angles. 

This paper is organized as follows. In section 2, we give a brief review 
on fixed points on the Zq-1 orbifold and the corresponding twisted states. 
Also, their selection rules for allowed Yukawa couplings and the strength of 
Yukawa couplings are reviewed. In section 3, we study systematically the 
possibility for realizing realistic lepton masses and mixing angles by using 
only renormalizable Yukawa couplings derived from Zq-1 orbifold models. 
We assume one pair of up and down sector Higgs fields. We consider the 
Dirac neutrino mass scenario in section 3.1. In section 3.2, we perform the 
same analysis for the seesaw scenario. Section 4 is devoted to conclusion and 
discussions. 


2 Orbifold models and selection rule 

Here we give a brief review on heterotic orbifold models, that is, the structure 
of fixed points on orbifolds, the selection rule for allowed Yukawa couplings 
and their Yukawa coupling strength. In particular, we concentrate our at¬ 
tention on Z 6 -1 orbifold models. (See for details of generic Z^ orbifolds 
Refs, m 21]-) 

2.1 Fixed points and twisted sectors 

An orbifold is defined as a division of a torus by a discrete rotation, i.e., 
a twist 6. The 6 D Z 6 -I orbifold is obtained by dividing T 6 by the twist 9, 
whose eigenvalues are diag(e 27r4//6 , e 2 ”/ 6 , e 2?n / 3 ), that is, a direct product of 
two 2D Zq orbifolds and a 2D Z 3 orbifold. The 2D Zq orbifold is obtained 
e.g. through dividing K 2 by the G 2 lattice and its automorphism, that is, 
the Coxeter element of G 2 algebra, which transforms the G 2 simple roots e\ 
and e 2 , 

Oe 1 —» — e 2 , $e 2 —> 3ei + 2e 2 , (1) 

i.e., the Z & twist. Similarly, we can obtain the 2D Z 3 orbifold by dividing 
R 2 by the SU{3) root lattice and its Coxeter element, which transforms the 
SU(3) simple roots e$ and ee, 

9e 5 —> e$, dee ~ 65 — e6, (2) 
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that is, the Z 3 rotation. 

There are two types of closed strings on an orbifold. One is the untwisted 
string, which is already closed on a torus before orbifold twisting. The other 
is the twisted string, which is a localized mode we are interested in. A twisted 
string has the following boundary condition, 

X\a = 2vr) = {9 k X) i {a = 0) + v\ (3) 

where the shift vector v l is on the torus lattice A. That is, the center of mass 
of a 9 k twisted string is localized at a fixed point /, which is defined as 

r = (<> k fy + v\ ( 4 ) 

The fixed point / is represented by the corresponding space group element 
( 9 k , v *). The fixed points, which differ by lattice vectors, correspond to equiv¬ 
alent fixed points. That implies that ( 9 k ,v l ) is equivalent to ( 9 k ,v l + (1 — 
9 k ) A). Hereafter the fA-twisted sector is denoted by T k . 

The 2D Z 3 orbifold has the following three fixed points under 9, 

Sza.i = (0.0), gg = (2/3,1/3), <*£>,( 1/3,2/3), (5) 

in the SU (3) simple root basis, and each of these is represented by the space 
group element as 

9zl i : (M el )> (6) 

where n = 0,1,2, up to (1 — 9)Asu{3)- The corresponding twisted ground 
states are denoted by r/^ 1 ) with n — 0 , 1 , 2 . 

Similarly, we can obtain fixed points on the 2D Z§ orbifold. The 9 twisted 
sector on the 2D orbifold has only one fixed point, 

»£’., = (0.0), P) 

in the G 2 simple root basis. The 9 2 twisted sector has three fixed points, 

92,2 = (0,0), jg 2 = (0,1/3), 92,2 = (0,2/3). (8) 

Note that these fixed points are defined up to the G 2 lattice. For example, one 
can use the fixed point g^ 2 = (1,2/3), which is equivalent to (0,2/3). The 

corresponding three twisted ground states are denoted by \g 'z f . 2 )- However, 
not all of the three points 2 are fixed points of 9. While r/^ 2 is also a fixed 
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point of the twist 0, the other two fixed points g^% 2 and 2 are transformed 
to each other by 9. Since physical states are constructed as 0-eigenstates, we 
take linear combinations of states corresponding to 2 and 2 as j2El 2D] 



with the eigenvalues 7 = ± 1 , while the state 7 (^ 2 ) corresponding to the 
fixed point 2 is by itself a 0-eigenstate. 

The 0 3 twisted sector has four fixed points, 

9S,3 = (°.°). = (0.1/2), 

sg >3 = ( 1 / 2 , 0 ), 9 |>3 = ( 1 / 2 , 1 / 2 ). ( 10 ) 

Recall that these fixed points are defined up to the G 2 lattice. For instance, 
the fixed point g ^ 3 = (1,1/2) is equivalent to (0,1/2). Not all of the four 
points are fixed points of the twist 0. The 0-eigenstates for each G 2 part are 
obtained as 



where 7 = l,u/,u/ 2 with u = e 2m / 3 . 

A fixed point on the 6D Zq-1 orbifold is obtained as a direct product of 
three fixed points, each coming from one of the two 2D Z e orbifolds and the 
2D Z 3 orbifold. The corresponding twisted ground state is obtained in the 
same manner. The 0 twisted sector has the following ground states, 

1 9ze,i) ® Ifl4e,i) ® I^Z3 ,i)j (12) 

for i — 0,1, 2. The 0 2 twisted sector has the following ground states, 

10^6,2) ® 1^0,2) ® 10^3,2)) 

1^8,25 7> ® 1 9zlz) ® 1 ^ 3 , 2 ) > 

\9z1 2 ) ® \9zlv i) ® WzL): ( 13 ) 

\9zl^ 7) ® 7') ® 102,2)> 
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for 7 , 7 ' = ±1 and j — 0,1,2. The 9 3 twisted sector has the following ground 
states, 

I^Z6,3) ® Wzl,3)i 
\9zl&l) ® I^S,3)> 

I^s) ® \9zlv'y')> ( 14 ) 

1^2,35 7) ® 1 9zl,^i) t 

where 7 , 7 ' = l,u;,u; 2 . 

2.2 Selection rule 

Here we give a brief review on the selection rule for Yukawa couplings in 
orbifold models. (See for their details Refs. [20, EJ.) The fixed point / of 
the 9 k twisted sector is denoted by its space group element, {9 k , (1 — 9 k )f), as 
said in the previous subsection. Thus, the three states corresponding to the 
three fixed points ( 6 ki , (1 — 9 ki )fi ) for i — 1, 2, 3 can couple if the product of 
their space group elements Y\ t {9 kl , (1 — 9 ki )fj) is equivalent to identity. That 
implies the space group selection rule for allowed Yukawa couplings requires 

m 

(9 ki , (1 — 9 ki )(fi + A)) = ( 1 , 0 ), (15) 

i 

because the fixed point ( 9 k , (1 — 9 k )f ) is equivalent to ( 9 k , (1 — 9 k )(f + A)). 
This space group selection rule includes the point group selection rule, that 
is, the product of twists must be identity, IL 9 ki = 1. In Z 6 -I orbifold models, 
the point group selection rule and //-momentum conservation m allow only 
the following couplings [25] , 

T^f;, t 2 f 2 T 2 . (16) 

The space group selection rule for the 2D Z 3 orbifold is simple. For a 
T 2 T 2 T 2 coupling, three states corresponding to fixed points l 2 , 9z% an<4 

9zl\, can coll pl e when the following equation is satisfied, 

ii + h + h — O (mod 3). (17) 

Also, for a TiT 2 T 3 coupling, three states can couple when the fixed point of 
Tf is the same as that of T 2 . Thus, the space group selection rule of the 2D 
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Z 3 part allows only diagonal couplings, that is, this part is not relevant to 
our purpose of deriving realistic Yukawa matrices with non-vanishing mixing 
angles in the case with the minimal number of Higgs fields. For a while, we 
will assume that all relevant states correspond to the same fixed point on the 
2D Z 3 orbifold. There is another possibility for a T 2 T 2 T 2 coupling, that is, 
the three states correspond to different fixed points on the 2D Z 3 orbifold, 
and in this case, we have a suppressed Yukawa coupling depending on the 
volume of the 2D Z 3 orbifold. We will give a comment on it later. 

We can discuss the space group selection rule for the 2D Z§ part. For 
T 2 T 2 T 2 couplings, the space group selection rule on the 2D Z e orbifold is 
exactly the same as that on the 2D Z 3 orbifold, i.e. eq. (ED, when we consider 
the basis of twisted states corresponding directly to fixed points. However, 
in Zq-1 orbifold models, we take linear combinations as in Eq. ® . The space 
group selection rule for T{T 2 T^ couplings is non-trivial. All couplings on 
the 2D Zq orbifold are allowed by the space group selection rule, because 
(1 — 9 )A = A. Thus, off-diagonal couplings are allowed for T 1 T 2 T 3 couplings 
in the 2D Z 6 orbifold part. Furthermore, the selection rule requires that the 
product of ^-eigenvalues of the coupling states must be equal to identity, that 
is, n 7 = 1. Therefore, the twisted states, which are relevant to our purpose, 
are the single Ti state, 



(18) 


and the Eve T 2 states, 



(19) 


where 7 = ±1, and the six T 3 states, 



( 20 ) 
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where 7 = l,u,u 2 . We have omitted the index for the fixed point on the 2D 
Z 3 orbifold, because we assume all states sit on the same fixed point on the 
2D Z 3 orbifold for the moment. 


2.3 Yukawa couplings 

The strength of a Yukawa coupling has been calculated by use of 2D confor¬ 
mal held theory. It depends on distances between fixed points. The Yukawa 
coupling strength of the T\T 2 T 3 coupling in Z 6 -1 orbifold models is obtained 
for the G2 x G2 part as 0001211 


Y = exp[-—f^ 3 Mf 23 ], 

/23=/2 — /3+A 

up to an overall normalization factor, where 

-§ R\ 0 0 \ 

3 R\ 0 0 

n p2 3 p2 ) 

U it2 2 **'2 

0 -§ R\ 3 R\ j 


l R\ 

3 p2 

2 W 

0 

0 


( 21 ) 


( 22 ) 


in the G 2 x G 2 root basis. Here, / 2 and f 3 denote hxed points of the T 2 and 
T 3 sectors, respectively, and Rj corresponds to the radius of the i-th torus, 
which can be written as the real part of the i-th Kahler modulus T t up to a 
constant factor. Here we follow Ref. EB for the normalization of i?*. (See 
Ref. [2E3 for the normalization of the moduli such that the transformation 
Ti —> Ti + i is a symmetry.) The imaginary parts of T % also contribute to mass 
matrices, i.e., eigenvalues, mixing angles and CP violating phases. However, 
here we consider only the real parts R L for simplicity. Since the states with 
hxed points in the same conjugacy class contribute to the Yukawa coupling, 
we take summation of these contributions in eq. m- However, the states 
corresponding to the nearest hxed points (/ 2 , f 3 ) contribute dominantly to 
the Yukawa coupling for a large value of R t . Hence, we calculate Yukawa 
couplings only by the contribution due to the nearest hxed points (/ 2 ,/s). 
Indeed such an approximation is valid when Ri is sufficiently large as in the 
cases we will study in the following sections. 

Similarly, the strength of T 2 T 2 T 2 Yukawa couplings is obtained in the 




( 23 ) 


basis of twisted states corresponding directly to fixed points as 
V = exp[-^/ 2 r 3 M/ 2 3], 

/23=/2 — /3+A 

where M is the same matrix as Eq. Here, / 2 and denote two of the 
three hxed points in the T 2 sector. Recall that when we choose two states, 
the other state, which is allowed to couple, is uniquely hxed in the basis of 
states corresponding directly to hxed points. 

Here, we give a comment on the Kahler metric. For a given fc, the 
6 k twisted states have the same Kahler metric, even if they correspond 
to different hxed points. Thus, the Kahler metric is irrelevant to mass 
ratios or mixing angles when the three families of leptons with the same 
577(3) x SU(2) x U (l)y quantum numbers are assigned to states in a single 
Tk sector. Indeed this type of assignment is required by the point group 
selection rule in order that non-vanishing mixing angles can be realized. 

3 Lepton masses and mixing angles 

Here we systematically study the possibility for deriving realistic lepton 
masses and mixing angles by use of twisted states, their selection rules and the 
Yukawa couplings, which are shown in the previous section. We assume that 
we obtain the standard gauge group SU{3 ) x SU(2 ) x U(1)y and three fami¬ 
lies of leptons in a Z e -l orbifold model. Indeed it is quite a nontrivial issue to 
construct a realistic model including the gauge group SU(3 ) x SU{2 ) x U (l)y 
and three families of quarks and leptons, but here we just assume them, be¬ 
cause our purpose is not to construct an explicit model, but to show the 
possibility for leading to realistic mixing angles and mass ratios. We also 
assume one pair of up and down Higgs fields. 

Charged lepton masses are well known and shown in Table |T| ! 23 • The 
table also shows the neutrino mass squared differences and mixing angles, 
which are consistent with recent experiments on neutrino oscillations |Tj. 

3.1 Dirac neutrino mass scenario 

First we consider the Dirac neutrino mass scenario. The relevant terms in 
the superpotential are 

W D HMYjijNj - H.L^YejijGj, (24) 
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Parameter 

Experimental value 

m e 

(0.51099892 ± 0.00000004) MeV 

"V 

(105.658369 ± 0.000009) MeV 

m T 

(1776.99 ± 0.275) MeV 

Arr4 

(8.1 ±0.3) x 10“ 5 eV 2 

Am§i 

(2.2 ±0.3) x 10“ 3 eV 2 

sin 2 9 12 

0.30 ±0.025 

sin 2 9 23 

0.50 ±0.065 

sin 2 0 13 

< 0.014 


Table 1: Experimental values of lepton masses and mixing angles. 


where H u and Hd are the up- and down-sector Higgs fields, L t are the SU(2) 
doublet leptons, W are the gauge singlet neutrinos, and are the SU(2) 
singlet charged leptons. After replacing the Higgs fields with their VEVs, v u 
and Vd, leptons gain the mass terms, 


W D —i'i(m D )i j Nj - e i (m e ) ij ej, (25) 

where (z^e*) = L i: ( m D )ij = ( v u Y v )ij , and (m e )p- = (v d Y e ) i:j . The mass 
matrices can be diagonalized as 

m D = m e = U e mf ag V}. (26) 

In terms of these unitary matrices, the lepton mixing matrix is defined as 

U M m = Uju:. (27) 

We consider all possible assignments of L i) N il e\ and H U)d to the twisted 
states shown in the previous section. The point group selection rule implies 
that we can obtain non-vanishing mixing angles only when each species of 
Li, Ni, and belongs to the same twisted sector. The Tf sector has a single 
state for the T 2 /Z 6 x T 2 /Z 6 part. Hence, we have to assign L % , Ni and 
to T 2 or T 3 to obtain non-vanishing mixing angles, while the Higgs fields can 
be assigned to T\. Therefore, there are five classes of assignments, which 
are shown in Table [21 For each of possible assignments, we examine the 
selection rule for allowed Yukawa couplings and calculate Yukawa couplings 
as functions of Ri and i? 2 - Then, varying these two parameters R\ and i? 2 , 
we try to fit the charged lepton mass ratios m e /m r and m^jmr, and the 
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Class 

L 

N 

e c 

H u 

H d 

Assignment 1 

T 2 

T 3 

T 3 

Ti 

T\ 

Assignment 2 

T 3 

f 2 

t 2 

T\ 

T\ 

Assignment 3 

t 2 

T 3 

t 2 

Ti 

t 2 

Assignment 4 

T 2 

T 2 

t 3 

To 

Ti 

Assignment 5 

T 2 

f 2 

T 2 

t 2 

T 2 


Table 2: Five classes of assignments 


ratio of neutrino mass squared difference and mixing angles $ 12 , 

$23, $ 13 - For the ease of presentation, we display the ratios of experimental 
values, which can be immediately obtained from Table |T| 


(rn e /m T ,m^/m T ) = (0.000288,0.0595), 
Am^/Am^ = 27. 


(28) 

(29) 


In particular, we are interested in deriving their orders, but not precise val¬ 
ues, because those values are obtained at the string scale. However, it is 
quite non-trivial to fit these six observables, m e /m T , m M /m r , Am^/Am^ 
sin 2 $i 2 , sin 2 $ 23 , sin 2 $13, only by two parameters, R\ and R 2 . It may be most 
important to realize the mass ratios of charged leptons m e /m T and m^/m T , 
because their values are measured very precisely by experiments. Hence, 
the two parameters, R\ and i? 2 , are almost fixed in order to fit m e /m T and 
m M /m T . In all the five classes of assignments, we can find cases that fit the 
charged lepton mass ratios, but not in all of them can we fit the neutrino 
oscillation data because mass matrix patterns are limited. 

We show our results in what follows. 


Assignment 1 

In this class of assignments, we can fit only the mass ratios m e /m T and 
m^/m T properly. For example, we take (R\,R\) = (26,33) in the following 
assignment, 


(L 1 ,L 2 ,L 3 ) = (f 2 (1) ,f 2 (2) ,f 2 (4 ’ 1) ), 


'li c 2> c 3. 


= (Tf.Tf ,f 3 ( 


(N 1 : N 2 ,N 3 ) = (fP.ff.ff), 

= (30) 


Then, we can realize the experimental values of m e /m T and m^/m T . How¬ 
ever, the neutrino oscillation data cannot be well accommodated. In this 
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sin 2 #12 = 5 x 10 5 , 


particular assignment, we obtain 


A ^31 

Am 2 21 


= 110 , 


sin 2 #23 = 6 x 10 3 , sin 2 # 13 = 6 x 10 


-li 


( 31 ) 


Thus, the mass squared difference ratio Am^/Am^i is large, and what is 
worse is that the mixing angles #12 and #23 are too small. This is the character 
of this class of assignments. Namely, we can fit the charged lepton mass 
ratios m e /m r and m^/m T , but the mixing angles are too small except the 
trivial cases. This character is shown in Table [3 Here and hereafter, by 
a trivial result we mean that at least one of charged leptons is massless, 
charged leptons or neutrinos have degenerate masses, or the mixing matrix 
is trivially the identity matrix. 


Assignment 2 

In this class of assignments, we can fit only the mass ratios m e /m r and 
m^/m T properly. For example, we take (_R 2 ,i? 2 ) = (21,26) in the following 
assignment, 


(L^La)^ 


_ ('T’(l) ^( 2 ) 'f 


(N u N. 1 ,N 3 ) = (fi l \fp,T^ 1) 


), 


(e;,eS,e§) = (T 2 « Tf.lf 1 ’), (H„H d ) = (T 1 ,T 1 ). 


(32) 


Then, we can realize the proper mass ratios m e /m T and The neutrino 

masses and mixing angles are not so satisfactory, however. We obtain 

sin 2 #12 = 4 x 10” 7 , 

sin 2 #13 = 9 x lO” 8 . (33) 

Thus, the mass squared difference ratio Am^/Am^ is large, and the mix¬ 
ing angles #12 and #23 are too small. That is the character of this class of 
assignments. Namely, we can fit the charged lepton mass ratios m e /m T and 
m /t /m r , but the mixing angles are too small except the trivial cases. This 
character is shown in Table 01 


Am h 

Am|i 


= 150, 


sin 2 #23 = 2 x 10 


-3 


Assignment 3 

In this class of assignments, we can fit only m e /m T , m M /m T , and Am^/A 777 J 1 
properly. For example, we take (i? 2 ,i?|) = (205,481) in the following assign¬ 
ment, 


(Li, L 2 , T 3 ) 


(rp(2) rp( 3) rr)(4,l)N 
\ 1 2 > J 2 )i 


(N h N 2 ,N 3 ) 


/■'T'CO 7 ^( 2 ) ffi(4,l)\ 
l J 3 5 J 3 1 1 3 )i 


12 




(e;, el, e|) = (if, if\ if 1 *), (H„ H d ) = (T„ If') 


(2b 


(34) 


Then, we can realize the proper mass ratios m e /m T and m tl /m T . On the 
other hand, we obtain 


= 

sin 2 003 = 


19, 

1 . 0 , 


sin 2 0 12 = 3 x 10 4 , 
sin 2 0 13 = 3 x 10~ 10 . 


(35) 


Thus, the mass squared difference ratio is compatible with the 

data, but the mixing angle 0 i2 (023) is too small (too big). The sizes of 
sin 2 012 , 23,13 vary from case to case, but they turn out to be either too big 
> 0.9 or too small < 0.1 in most cases. Sometimes they are 0.3 or 0.5, 
but no case leads to the right values of the three mixing angles at the same 
time. That is the character of this class of assignments. Namely, we can fit 
the charged lepton mass ratios m e /m T and and the neutrino mass 

squared difference ratio, but the mixing angles are either too 

big or too small. This character is shown in Table 01 

Assignment 4 

In this class of assignments, we can fit the neutrino oscillation data as 
well as the mass ratios m e /m T and m^/m T . For example, we take (i? 2 , R\) = 
(26, 21 ) in the following assignment, 


(L i, L2, L 3 ) 


rp(3) rp{ 4 ,- 1 ) 


1 


{N u N 2 ,N 3 ) 


(t 2 ( 2 ) ,t 2 (4 ’ 1} 


— 1 )\ 

) 1 2 )i 


( e li e 2i e V) ~ 




), 


(H u ,H d ) = (f?\f 1 ). 


(36) 


Then, we can realize the experimental values of m e /m T and m /t /m T . More¬ 
over, in this assignment we obtain 


A ml 

- — = 14 

Am 2 ’ 


21 


sin 2 0 23 = 0.70, 


sin 2 0 i 2 = 0.38, 
sin 2 0 i 3 = 6.3 x 10~ 6 . 


(37) 


Thus, this assignment can realize orders of six observables only by two param¬ 
eters. Almost the same result can be derived from the following assignment, 


(Li,L 2 ,L 3 ) = (T 2 2) , T 2 3) , T 2 4 ’ _1) ), (N h N 2 ,N 3 ) = {T?>,T?' L) ,T, 


( 2 ) rp(4,l) rp{ 4 ,- 1 ) 


_ /T 1 ! 2 ) TA) TM’-O 


(e?,e£,e§) = 


n"“), (H u ,H d ) = (tp,T 1 ). 


). 

(38) 
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There are also several other assignments leading to similar results with smaller 
values of Am^/Am|i in the range of 


A nn\ x 


< 14. 


In Table 0 the best fitting result is shown. 


(39) 


Assignment 5 

In this class of assignments, we can fit only m e /m T , m^/m T , and 
properly. For example, we take (Rf, R 2 ) = (243,419) in the following assign¬ 
ment, 


(L i, L 2 , L 3 ) 

( e l) e 2, e 3) = 


( r v0-) r p(2) 7 k 4 ) \ 
l J 2 J i 2 J i 2 b 



44 3) 
? 4 2 


j 4 2 b 


(N U N 2 ,N 3 ) = (f«, if, if 1 ’), 
(H u ,H d ) = (fP,ft' ) ). (40) 


Then, we can realize the proper mass ratios m e /m T and m^/m T . On the 
other hand, we obtain 


A/noi 

a^4 =28 ’ 

sin 2 6*23 = 0.06, 


sin 2 6*12 = 1 x 10 3 , 
sin 2 6*13 = 3 x 10 -14 . 


(41) 


Thus, the mass squared difference ratio Am 31 / Amh is compatible with the 
data, but the mixing angles 6*12 and 6*23 are too small. The size of sin 2 6*23 
varies from case to case, but it turns out to be either too big ~ 0.9 or too 
small ~ 0.1 in any case. That is the character of this class of assignments. 
Namely, we can fit the charged lepton mass ratios m e /m T and m M /m T , but 
$12 is too small and $23 is either too big or too small. This character is shown 
in Table 01 


As a result, we can realize the charged lepton mass ratios, the ratio of 
neutrino mass squared differences and the lepton mixing angles in certain 
cases of Assignment 4, but not in the other classes of assignments. Let us also 
note that the neutrino mass spectrum has normal hierarchy with vanishing 
or very small lightest neutrino mass in every case leading to reasonable value 
of Am^/Amli- 

So far, we have assumed that all states correspond to the same fixed point 
on the 2D Z 3 orbifold. However, that case leads to ( 1^)33 = 0(1) which is not 
realistic. We need a suppression factor of O(10 -12 — 10~ 13 ) to fit the overall 
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Class 

Am 2 31 /Am 2 21 

sin 2 0 i 2 

sin 2 023 

sin 2 0i 3 

Assignment 1 

~ 100 

<10" 5 

<io - 2 

< 10-' 

Assignment 2 

- 100 

<10' 5 

<io - 2 

<10' 7 

Assignment 3 

>1-4 




Assignment 4 

14 

0.38 

0.70 

6.3 x 10" 6 

Assignment 5 

~ 28 

< 0.09 


<10' 2 

Experimental values 

27 

0.30 

0.50 

0.000 


Table 3: Characteristics of each assignment 1-5 in the Dirac neutrino case. 
Typical behavior of each value is described for combinations resulting in 
relatively good fits in a given assignment except Assignment 4. The row 
corresponding to Assignment 4 shows the best fit. We omit m e /m T and 
m /t /m T because they can be fit in all the assignments. 


magnitude. In Assignment 4, the neutrino Yukawa couplings originate from 
T 2 T 2 T 2 couplings. On the 2D Z 3 orbifold, Yukawa couplings corresponding to 
three different fixed points are allowed, too. Here, we assume that Lj, N t and 
H u sit at three different fixed points on the 2D Z 3 orbifold. In this case, the 
neutrino Yukawa couplings universally have an exponential suppression factor 
like eqs.(J2U), (1221) and 423). Then, we can derive (W) 33 = 0(10 12 ), when 
we take the radius of the third torus as R\ = 0(1000). That implies that the 
compactihcation scale is smaller by O(10 _1 —10” 2 ) than the string scale. (See 
for phenomenological aspects of such a scenario e.g. Refs. uanzmang.) 

3.2 Seesaw scenario 

In this subsection, considering the seesaw scenario of the neutrino mass ma¬ 
trix, we systematically carry out the same analysis as in the previous sub¬ 
section. The relevant terms in the superpotential are 

W D H u L i (Yl)i :j N j - (42) 

where the third term is the right-handed majorana neutrino mass term. After 
replacing the Higgs fields with their VEVs, we obtain the mass terms, 

W D - e i (m e ) ij e c j - (43) 
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We assume M N ^ rrio 1 and in this case the heavy fields N t can be integrated 
out. Then the lighter degrees of freedom have the effective mass terms, 



(44) 


where 


rn u 


— 7TIdM n 1 7TI 


T 

O' 


(45) 


In string models, the natural order of mass is the string scale. On the 
other hand, the mass scale of M^ is phenomenologically required to be an 
intermediate scale between the string scale and the weak scale. Such an 
intermediate scale may be obtained by VEVs of some fields, as a comment 
will be given at the end of this section. However, such a scenario is quite 
model-dependent. Thus, we assume that M N is proportional to the identity 
for simplicity. Under such an assumption, we obtain 


M N l m D m T D = , 


(46) 

(47) 


diag 

V 


where U u is the same diagonalizing matrix as the one in the Dirac neutrino 
mass scenario. Hence, the mixing angles are almost the same as those in the 
Dirac neutrino mass scenario. On the other hand, neutrino mass hierarchy 
is enhanced compared with that in the Dirac neutrino mass scenario. The 
results are as follows. 

Assignment 1 

Like the Dirac neutrino mass scenario, we can fit only the charged lepton 
mass ratios properly. The neutrino mass squared difference ratio is too large, 
and mixing angles are too small. Those results are shown in Table 0], 

Assignment 2 

We can fit only the charged lepton mass ratios. The neutrino mass 
squared difference ratio is too large, and mixing angles are too small, as 
shown in Table 0] 

Assignment 3 

We can fit the charged lepton mass ratios. However, the neutrino mass 
squared difference ratio turns out to be rather large > 50 except for the case 
that gives = 2. Also, the mixing angles are typically either too 

big or too small. 
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Class 

Am 2 31 /Am 2 2l 

sin" 1 0i2 

sin 2 023 

sin 2 0 i 3 

Assignment 1 

~ 6000 

<10" 5 

<io- 2 

< 10 -' 

Assignment 2 
Assignment 3 

~ 7000 
> 2 

<10' 5 

<io- 2 

<10' 7 

Assignment 4 

29 

0.32 

0.48 

3.6 x 10" 6 

Assignment 5 

~ 28 

< 0.09 


<io- 2 

Experimental values 

27 

0.30 

0.50 

0.000 


Tabic 4: Characteristics of each assignment 1-5 in the seesaw case. Typical 
behavior of each parameter is described for combinations resulting in rela¬ 
tively good fits in a given assignment. The row corresponding to Assignment 
4 shows the best fit. We omit m e /m T and m M /m T because they can be fit in 
all the assignments. 


Assignment 4 

Like the Dirac neutrino mass scenario, we can properly fit both the 
charged lepton mass ratios and the neutrino oscillation data. For example 
we take (i?f, K\) = (23, 26) in the following assignment, 


(L i, L2, L 3 ) 


2^(2) ^(( 4 ,1) 


(N u N 2 ,N 3 ) 


/^( 2 ) 3 ) ^( 4 , 1 ) 


(e c v e 2,e c 3 ) 


^( 2 ) rj t ( 4 , 1 ) 


), 


(H u ,H d ) = (T?\f 1 ). 


(48) 


Then, we can realize the experimental values of m e /m T and m M /m T . Fur¬ 
thermore, in this assignment we obtain 


A "4 

Amli 


= 29, 


Sill 


e 23 = 0.48, 


sin 2 9 12 = 0.32, 
sin 2 di 3 = 3.6 x 10 -6 . 


(49) 


Thus, this assignment can realize orders of six observables only by two pa¬ 
rameters. This class of assignments include several other cases leading to 
almost the same results. 

Assignment 5 

We can fit the charged lepton mass ratios and However, 0 12 

is too small and 6 23 is either too big or too small, as shown in Table 01 

As a result, we can realize the charged lepton mass ratios, the ratio of 
neutrino mass squared differences and the lepton mixing angles in certain 


17 






cases of Assignment 4, when we assume that the right-handed majorana 
neutrino mass matrix Mjv is proportional to the identity matrix. As was in 
the Dirac scenario, the neutrino mass spectrum has normal hierarchy with 
vanishing lightest neutrino mass. The neutrino Yukawa coupling (Y v ) 33 is 
0(1). We need Mn = O(10 15 ) GeV for the seesaw mechanism to produce 
the correct neutrino mass scale. 

I 11 the above analysis, we have assumed (M/v)^ = for simplicity. It 

is nontrivial to realize such a right-handed neutrino majorana mass matrix. 
Majorana mass terms can be generated by the (n + 2)-point couplings 

YijNiNj^ ■ ■ ■ 0 n ), (50) 

after 0j develop their VEVs.' Thus, the right-handed neutrino mass matrix 
depends on VEVs of scalar fields and the selection rule for higher dimensional 
operators HQ, that is, quite a model-dependent feature. In certain cases, we 
may obtain (Mjv)jj = However, we, in general, obtain nontrivial ma¬ 

jorana matrix ( Mm)^ if it is derived from YijNiNj^cpi ■ ■ ■ (p n ) through VEVs 
of (0i • • • 4> n )- hr such a case, assignments other than assignment 4 might 
lead to realistic results. Thus, it would be interesting to study such case 
somehow systematically. However, that is beyond our scope. We leave it for 
future study. 

4 Conclusion 

We have systematically studied the possibility for realizing lepton masses 
and mixing angles by use of only renormalizable couplings derived from Z 6 -l 
heterotic orbifold models. We have assumed one pair of up and down type 
Higgs fields. We have found Assignment 4 has such a possibility in both 
the Dirac neutrino mass scenario and the simple seesaw scenario where the 
right-handed majorana mass matrix is proportional to identity. The resulting 
neutrino mass spectrum shows normal hierarchy. 

It is quite non-trivial to fit six observables only by two parameters R\ and 
i? 2 . However, how to stabilize these moduli at proper values is an important 
issue to study further. Moreover, if F-components of moduli fields contribute 
to SUSY breaking, models corresponding to Assignment 4 may show a certain 
pattern of soft SUSY breaking parameters. It is interesting to study their 
effects on flavor violation. (See e.g. [22] and references therein.) 

7 Such VEVs may be given e.g. through anomalous U(l) breaking m- 
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In addition to the real parts Ri of the moduli 7), mass matrices also 
depend on imaginary parts of T t , although we have fixed them to vanish in 
our analysis. When such imaginary parameters are included, other types of 
assignments may lead to realistic results. Thus, it is interesting to carry out 
the same analysis including imaginary parts of T*. 

We have systematically studied all possible assignments of leptons and 
Higgs fields to twisted states, but not constructed an explicit heterotic orb- 
ifold model by fixing gauge shifts and Wilson lines. It is also important to 
construct explicitly heterotic models corresponding to Assignment 4. 

In addition, it is interesting to extend our analysis to other Zjy and Zjy x 
Z M orbifold models including the quark sector. In principle, systematical 
studies as in this paper are possible for other orbifold models. Such a study 
will be done elsewhere. 

Although we have discussed the possibility for obtaining realistic Yukawa 
matrices from only stringy renormalizable couplings in this paper, another 
possibility for realistic Yukawa matrices is that higher dimensional opera¬ 
tors play roles to generate effective Yukawa couplings after symmetry break¬ 
ing. (See e.g. P3J.) Particle mixing through symmetry breaking is also one 
possibility [33] • If would also be interesting to study such a possibility some¬ 
how systematically. In the former case, it is quite important to study sym¬ 
metries to control higher dimensional operators in string models. 
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